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Abstract 

A large class of first order partial nonlinear differential equations 
in two independent variables which possess an infinite set of polyno- 
mial conservation laws derived from an explicit generating function 
is constructed. The conserved charge densities are all homogeneous 
polynomials in the unknown functions which satisfy the differential 
equations in question. The simplest member of the class of equa- 
tions is related to the Born-Infeld equation in two dimensions. It is 
observed that some members of this class possess identical charge den- 
sities. This enables the construction of a set of multivariable equations 
with an infinite number of conservation laws. 







1 Introduction and examples. 



The purpose of this article is to demonstrate the existence of a class of non- 
linear first order differential equations in 2 independent variables, but an 
arbitrary number of dependent ones which admit an infinite number of in- 
dependent polynomial conservation laws. There are general results which 
guarantee the existence of an infinite set of conservation laws for evolution 
equations in two dimensions but there are comparatively few recorded 
cases where the conserved densities may be written as explicit polynomials 
in the dependent variables. One purpose of this article is to demonstrate a 
class of equations in which these conservation laws may be derived from the 
use of a generating function, another is to exploit these results to exhibit an 
example with polynomial conservation laws in more than two independent 
variables. 

The equations considered are in fact particular examples of the equa- 
tions of hydrodynamic type studied by the Russian School of Novikov and 
Dubrovin 000], but the exphcit nature of the construction of the conserva- 
tion laws presented here may be useful. For a flavour of the results, consider 
first the following examples for 4 unknown functions Ui, i = 1 ... 4 which 
were discovered by a combination of guesswork and computer verification. 
These results led to a generalisation which will be described in two stages, 
first using symmetric polynomials, then without symmetry. General proofs 
will be given in the final two sections. 



1.1 Example 1 

Consider the following first order differential equations 

dui dui 



du2 , .dU2 , . 

— = (n, + ..4 + ..3)^ (1.1) 

duj, du3 
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Construct the matrices M[n,p, g,r] as follows: 



M[n,p, q, r] 



U? Un Ua 



p p 



Ml Uo 



P 



p 



U\ Un Un U\ 



[1.2) 



and define A =^ det M[3, 2, 1, 0]. 

Then the conservation laws for ( |1 . 1| ) are simply 

d /detM[n,2,l,0]\ d /det M[n, 3, 1, 0] \ 



dt 



A 



dx 



A 



V positive n > 3. (1.3) 



The conserved densities are clearly linearly independent. 
1.2 Example 2. 

Consider instead the following first order differential equations given by 

dui , , du\ 

{U4U2 + U2U3 + U^Uij- 



dt 



dx 



:i.4) 



together with its cyclic replacements. Then the conservation laws for this 
system of equations are simply 



d /detM[n,2,l,0]\ _ d /detM[ri,3,2,0]\ 
dt [ A ) [ A J ■ 

1.3 Example 3. 

Change once again the first order differential equations as follows; 

dui dui 



:i.5) 



U2USU4- 



dt dx 
Now the conservation laws for flOD become 



together with cyclic replacements. 



:i.6) 



d /detM[n,2,l,0]~ 
dt \ A 



d /detM[n,3,2,l]~ 
'dx \ A 



:i.7) 
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The general pattern becomes clear. The conserved charge density is a 
symmetric polynomial. It is the same one for each example. This leads to 
the evident consequence that a correponding result holds good for a linear 
combination of equations (|1.1|) , ( |1.4| ) and (|1.6| ). Furthermore, one may take 
equation (|1 . 1|) with independent variables x, t and add to the right hand 
side that of ( [L.4| ) with independent variables y, t and ( |1.6| ) with independent 
variables z, t and thus construct an example of a set of partial differen- 
tial equations in 4 dimensions, with independent variables x, y, z, t which 
admits an infinite number of conservation laws. This construction clearly 
generalises to any number of variables. This point will be elaborated later. 



1.4 Implicit Solutions. 



Example ( |1.1| ) can be integrated implicitly in the case of two dependent 
variables. In the case of three functions ui, U2, U3, the equations ( |1.4] ) 
can be reduced by the hodograph method of interchanging the roles of the 
independent variables x, t and the dependent variables Mi, U2 to give 

dx dt 
-U2U3- 



du2 du2 

dx dt 
-U3U1- 



dui dui 
dus dt Ui{u2 — Us) du3 dt 



dui du2 U2{us — ui) du2 dui 

Introduce three arbitrary functions of a single variable, /("Ui), g{u2), and 
F{z) where 

u'f{u')du'+ / v'g{v')dv'. (1.9) 



Then the solution of (|1.8| ) is given explicitly by quadratures as 

«l rU2 r T'^ u'f{u')du'+ v'g{v')dv' 

f{u')du'+ g{v')dv' - F{z)-^dz 



/Ml /•«2 r f"^ u' f{u')du'+ v'g(v')dv' 

u'^f{u')du'+ v'^g{v')dv' - F{z)dz 

U3 = F{ u'f{u')du'+ v'g{v')dv'). (1.10) 
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2 Application. 



In the case of two dependent variables, ui, U2 equations of the form 
arose in a study of the Born-Infeld equation with J.MulveyP]. The argu- 
ment may be reversed as follows; Suppose a second order Lorentz invariant 
equation (regarding t, x as light-cone co-ordinates) is sought by setting 

„, _ J\dtdx) 

«2 = (2.1) 

Then the resulting equations to determine /, g are so restrictive that essen- 
tially the only solution reproduces the Born-Infeld equation, which in 
dimensions is given in terms of the parameter A by 

^dx' dt^ ^dt' dx^ ^ dxdt'dxdt ■ ^ ' ' 

and Ml, U2 are roots of the characteristic equation 

3 First Generalisation 

Let Sn be the nth elementary symmetric polynomial of degree n (i.e. the one 
in which no variable appears more than linearly in any term) in variables 
Uk, k = 1 . . . N, n < N and S"^ be the nth elementary symmetric polynomial 
of degree n in iV — 1 variables Uk, k = l...,j,... N, where the superscript 
j denotes that the variable uj is omitted. Postulate the existence of the 
differential equations 

?^^s;^^ .^i,...,iv. (3.1, 

Introduce the functions (pk, k = . . . N given by 

^° = ni=f (A - u,) 
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and the recurrence relations 



(l)k + Mk-\ = —^nT^ ~T' (3-2) 



Sk 

with = 1. Then the result is that in consequence of the equations (|3^ 
the following relation holds 

<90o 



dt dx ■ ^"^-"^^ 

Comparison of the coefficients of inverse powers of A on both sides of a 
formal expansion of this equation gives the various conservation laws, since 
the quantities (j)^, (pn are generating functions for the conserved densities 
and currents respectively. These coefficients will be symmetric polynomials 
in the dependent variables which can in fact be identified with the evident 
extension of the ratios of the determinants of Vandermonde type presented 
in the three examples. The proof of these results hinges upon the fact that 
it is often easier to construct a generating function for a set of conservation 
laws than the conservation laws themselves H. It runs as follows; 



3.1 Proof 

The relations ( p.2|) can be solved in terms of 0o to yield 

r=k 

with 5*0 = 1. Then using the facts that 



(3.4) 



dS 

= "^Li; (3.5) 

^i'^n-l = Sn — S^, (3.6) 
which are easily established, the following relations may be deduced; 

rifh f) ^ o 

<^Vn / / V^/ 1 \n-r c \n-r , r=0 \'-'"'J 

a^-^'^^^a^.V^^^ + x^, 
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A — Uj dx 

^'^^ (3.7) 



dt 



Here relations (|3.5|) , (|3.5|) and equations have been used. This is just 

the result ( ^.3[ ) claimed. 

The connection with the specific examples written in determinantal form 

P , , . ^ .1 , . detM[n,2,l,0l 
requires some lurther explanation. Consider the quantity ^ . 

This is a symmetric polynomial of degree n — 3 in the variables ui, U2, us, U4 
such that the coefficient of every term is unity. Now consider the expansion 
in powers of of the expression 

^ {x-u,){x-u2){x-us){x-u,y ^^-^^ 

Then the coefficient of A"" in the power series expansion of this expression 

detMk,2,l,0] ^, 
IS the quantity ^ . This may be seen m various ways, the most 

direct of which is the observation that the coefficients in the expansion of each 
individual factor are all unity, and in the multiplication of the geometrical 
series all possible terms arise with unit coefficient. It can also be seen from 
the definition of the matrices involved, (|1.2| ) with the help of the expansion 
in partial fractions 

1 1 



(A - Ui){X - U2){X - %)(A - U4) (A - Ui){Ui - U2){Ui - U3){ui - U4) 

1 1 



(A - U2){U2 - Ui){u2 - U3){U2 - U4) (A - " Ui){u3 - U2){u-i - M4) 

^ {X - Ui){Ui - Ui){Ui - U2){U4 - U^i) ' ^^'^^ 

by considering the expansion of the determinant on the first row, which 
corresponds to taking the coefficient of A~" on the right hand side of this 

exp.e.io„. In a stoila. manner the identification of ^fli^MM ^i^l, 

the coefficient of A"" in the expansion of 

— A + Ml + ^2 + ^3 + ^4 



(A - Ui){X - U2){X - M3)(A - M4) 
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may be best accomplished by way of the identity ( p.9| ) and expansions of the 
determinant on the first row. The results for the other determinants can be 
established likewise. 



3.2 Multivariable case 

Since the same function 0o enters into the generating equation ( p.3|) for each 
n, these relations may be combined linearly to produce a multivariable ex- 
tension; Suppose there are independent variables t, Xi, i = 1 . . . N — 1 and 
the N dependent variables Uj, j = 1 . . . N satisfy the quasi-linear system of 
equations 

(3.10) 

Then this set of equations admits an infinite number of linearly independent 
polynomial conservation laws, obtained by equating the coefficients of inverse 
powers of A in a formal expansion of the generating function 

It is not clear at this stage whether the existence of these laws is sufficient 
to guarantee integrability of (|3.1CI|) . 



4 Further asymmetric extension 

The system admits a further generalisation which breaks the symmetry. Con- 
sider the following system of first order differential equations; 

dui dui 

— = {{ni - l)Ui + n2U2 + '^3U3)-Q^ 

— = (niui (n2 - 1)^2 + nsMs)-^ (4.1) 

^ = (niMi + ^2^2 + (ria - 1)^X3)^. 

This system possesses an infinite number of conservation laws each individual 
member of which corresponds to equating coefficients of A~'^ on both sides 
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of the equation 

d 1 d -\ + riiUi + n2U2 + ^3^3 



dt (A - Mi)"i(A - M2)"2(A - M3)"3 dx (A - Mi)"i(A - W2)"2(A - Ms)"^' 

(4.2) 

This example may be further generahsed by scahng ui, U2, to make the 
'diagonal' entries on the right of equations (|4.1|) arbitrary. If a solution 
of these equations of the form Ui = xfi{t), then these equations assume 
the typical form of dynamical evolution equations. The following general 
result may be abstracted from i ^A\} and the previous examples. Consider 
the expression 

M „ k=N J 

*M = '^N • (4-3) 

k=l 

The construction in the numerator is just a device to take the sum of only 

k=N 

the first M + 1 terms in the expansion of (A — Uk)"^^ in decreasing powers 

fc=i 

of A. It is assumed that if any of the Uj are not positive integers then the 
contour integral is taken round a contour including the origin, but excluding 
any other possible singularities. Then 



^ (A - «,)^(V ifliX 



dui Xi^'^HX-u^ \ dui P'qJ ^J-^ 2mz^+'^' 

^ I M k=N r 



AE%(A-^.,) t=Y 271*^^+' 



n ^ f / A(l-.)m=f(A-z..)"^ ^ 



'X^^^{X-Ui)J 2m{X-zUi) ^ ^ 

^-A/_o ^ function independent of A. 

dUi 



Note that the index i is free in (|4.4|) . In the last two steps first the geometric 
series is summed, then the contour integral performed yielding a function 
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which after removing a factor A *^ is independent of A. This function is, 
as can be seen from ( [4.4|) , equivalent to setting A = in the expression 

M k=N J 

These considerations give the final generalisation; The function $o serves 
as a generating function for the conservation laws for a very wide class of 
nonlinear differential equations as a consequence of the result 



A 



M 



M 



dt dx 

where the dependent functions Ui satisfy the equations 



(4.5) 



duj 
dt 



duj 
dx 



M 



k=N 



\=Ui 



(4.6) 



The individual conservation laws are obtained by equating powers of — on 

A 



both sides of ( [4.5|) . In the case where all ni are positive integers, this final 
generalisation is equivalent to the symmetric case when subsets of the vari- 
ables Ui are identified. If a solution with power law x dependence of the 
equations of the type ( ^^ ) is sought, then these equations become nonlinear 
evolution equations of the type studied in dynamical systems theory. For 
example, if Uj{x,t) = xvjif) then the x dependence of the equations ( |4.1| ) 
drops out. 



5 Coda 

A recent attempt to prove the conservation laws for example 1 directly led 
to the further generalisation; 




This equation is a conservation law for equations (|1 . 1| ) for any arbitrary differ- 
entiable functions /, h, k. This result may be proved by straightforward 
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differentiation after expansion on tlie first row of tlie determinants involved. 
Clearly this result will generalise. 
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